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Abstract. We provide a new dynamic approach to scenario generation for 
the purposes of risk management in the banking industry. We connect ideas 
from conventional techniques — like historical and Monte Carlo simulation — 
and we come up with a hybrid method that shares the advantages of standard 
procedures but eliminates several of their drawbacks. Instead of considering 
the static problem of constructing one or ten day ahead distributions for vectors 
of risk factors, we embed the problem into a dynamic framework, where any 
time horizon can be consistently simulated. Additionally, we use standard 
models from mathematical finance for each risk factor, whence bridging the 
worlds of trading and risk management. 

Our approach is based on stochastic differential equations (SDEs), like the 
HJM-equation or the Black-Scholes equation, governing the time evolution of 
risk factors, on an empirical calibration method to the market for the chosen 
SDEs, and on an Euler scheme (or high-order schemes) for the numerical eval- 
uation of the respective SDEs. The empirical calibration procedure presented 
in this paper can be seen as the SDE-counterpart of the so called Filtered 
Historical Simulation method; the behavior of volatility stems in our case out 
of the assumptions on the underlying SDEs. Furthermore, we arc able to 
easily incorporate "middle-size" and "large-size" events within our framework 
always making a precise distinction between the information obtained from 
the market and the one coming from the necessary a-priori intuition of the 
risk manager. 

Results of one concrete implementation are provided. 

Key Words: risk management, stochastic (partial) differential equation, 
calibration, historical simulation, time series, jump processes 



1. Introduction 

A core part of risk management in the banking industry is the identification of 
risk factors and the generation of sc;enarios for a one up to ten days time horizon. 
This task involves dealing with basically three important issues and technical re- 
quirements. First, the scenarios need to extract the market information from the 
available time scries of the risk factors. Mathematically speaking, this problem is 
tackled via the identification of random variables that describe the risk factors one 
or ten days ahead and that match the most important stylized facts in the time 
series. A major complication in this classical statistics (or econometrics) question is 
that the available time series are extremely short (about 500 business days) in com- 
parison with the number of risk factors (several thousands). Even if the time series 
were long enough, standard statistical methods would take an excessively heavy 
computational efi'ort to come up with reliable solutions for such a high dimensional 
problem. 

Second, the shortness of the time series implies that one has to separately model 
extreme events, or even middle size events, in order to complement the information 
provided by the time series and to allow for cautious prediction. The third issue 
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is that the procedure for scenarios generation should be computationally fast since 
they have to be generated every (business) day in order to be able to evaluate the 
risk to which all the portfolios of the institution are exposed. 

There are two families of ways to deal with the problem of scenario generation 
in risk management, namely the historical and the distributional approaches. Both 
strategies produce random variables from which one can sample the scenarios: 

Distributional or Parametric Approach. The scenarios are constructed by 
sampling a multivariate parametric distribution that has been calibrated using the 
historical data. One can also fit certain distributions to vectors of risk factors and 
then use a copula to put together all these lower dimensional models. There are 
two major difficulties associated to this method: first, most statistical estimators 
of the distribution parameters that need to be calibrated yield excessively wide 
confidence intervals for the estimated coefficients due to the relative shortness of 
the historical time series in comparison with the dimensionality of the problem. 
Second, the available historical data may not correspond to any of the standard 
frequently used parametric distributions. 

Historical Approach. The scenarios are obtained by directly sampling the his- 
torical data (eventually after rescaling) . This method is much privileged in practice 
due to its ease of implementation and because it does not require any assumption on 
the historical distributions of the risk factors. This does not mean that this method 
is hypothesis free; indeed, its statistical legitimacy comes from the so called boot- 
strap method (see for instance [1]), whose requirements are not always satisfied, 
the most important of them being the iid (independent and identically distributed) 
character of the risk factors returns. Additionally, this method is obviously inca- 
pable of predicting rare (extreme) phenomena that did not happen before in the 
past; consequently these rare phenomena need to be added "by hand" (which is 
not a disadvantage but simply means that the time series is too short to reflect all 
possible risks). Several strategies have been proposed to tackle these problems; a 
particularly sophisticated and powerful one is the so called filtered historical simula- 
tion [2| , where the historical returns are "filtered" by volatility functions calibrated 
to each time series using GARCH type processes. 

Empirical calibration of SDEs. The approach that we present in this article 
is reminiscent of both the historical and the distributional approaches and tries to 
"interpolate" between them. It allows us to properly capture the stylized facts on 
dependence of the given time series and to easily carry out additional distributional 
tuning. We create the random variable for future scenarios based on stochastic 
differential equations, which describe the local dynamics of the respective risk 
factors ( "distributional aspect" ) . The stochastic differential equations are standard 
models from mathematical finance like the Black-Merton-Scholes, Cox-IngersoU- 
Ross, Heath-Jarrow-Morton models, etc, which are adequately chosen for each in- 
dividual class of risk factors; in particular all these equations arc free of arbitrage. 
The use of risk management models which could also be used for pricing or hedging 
derivative contracts is a tremendous conceptual advantage. This way we bridge 
between the worlds of trading and hedging. In particular the approach allows to 
actually come up with conceptually coherent methods for risk management which 
could in principle also be used for trading purposes. 

We calibrate the SDEs not via an optimization process as it is customary, but by 
directly using the historical data to construct local characteristics that capture the 
stylized facts of the time series. The method relies on the central limit theorem for 
covariance processes. Under certain conditions that we specify later on in the 
paper the calibrated SPDEs converge properly if the observed data stem from 
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an SPDE (which is basically the law of large numbers behind the construction). 
Heavy tails are generated by distinguishing between "trading time" and "clock 
time" (randomization of the period of simulation) and by adding extreme events, 
possibly with heavy tailed distributions. The first effect reflects that the flow of 
information runs at different speeds under different market volatility conditions and 
is responsible for "middle-size" events; the second effect reflects the possibility of 
the occurrence of rare (extreme) events not present in the historical data. 

Having the filtered historical simulation approach (FHS) [2] in mind, the em- 
pirical calibration approach can be seen as the counterpart of FHS on the side of 
mathematical finance, i.e., instead of using the GARCH- methodology to capture 
volatility, we apply standard models from mathematical finance and random time 
changes for the same purpose (see Section [s] and |4|. To be precise, the issue of cap- 
turing volatility correctly is addressed by rescaling with respect to realized volatility 
parameters of the time series via a sliding window method and by imposing specific 
assumptions on the dependence of volatility on the underlying risk factors. 

2. Risk factors driven by diffusive stochastic differential equations 

In this section we outline the basic setting of Wiener driven stochastic (partial) 
differential equations, which underlie the time evolution of risk factors. For sto- 
chastic (partial) differential equations we refer to ^ as main reference. Details on 
term structure equations can be found in [B]. 

Our formal setup consists of a filtered probability space (fi, 5", (3^t)j>Q, P), where 
i^t)t>o is ^ complete filtration and P a probability measure. The basis is assumed 
to carry a d-dimensional standard Brownian motion B. We use as models for the 
risk factors stochastic differential equations from mathematical finance; this leads 
us to stochastic partial differential equations of the following general type: let H 
be the Hilbert space where risk factors take values, then we consider 

d 

(2.1) dYt - (/ii(Ft) + f^2{Yt))dt + J2 <^iyt) • KdBl, 

i=l 

(2.2) Yo e H, 
where 

(2.3) <y{Y) -.Ho^Ho 

is an invertible linear map on some closed subspace Hq C H containing the real 



4.1 



span (Ai, . . . , Xd) of the set of return directions Ai, . . . , Ad G -ff (see Remark 
for the use of the word return here). Such volatility structures are usually called 
constant direction volatilities, see for instance [S]. The "volatility factor" a is 
chosen appropriately for the respective risk factors in order to exclude immediate 
arbitrages (like, e.g., negative interest rates). It should be interpreted as an a priori 
given factor governing the shape of the support of the risk factors; it is hence a 
geometric factor. The vector field fii corresponds to appropriately chosen no- 
arbitrage conditions; fi2 is another vector field that lies in the span (Ai, . . . , Xd) and 
corresponds to an appropriate Girsanov type change of measure. 

We suppose that the usual conditions on the drift ni hold, that is, it can be 
written as 

(2.4) fi,{Y) = AY + ^I3{Y) 

for some Y G dom(j^), where A denotes the generator of a strongly continuous 
semigroup which in our case will always be a shift to the right or the identity 
semigroup on the Hilbert space of risk factors H. The maps fj.^, ^2 ■ H H are 
smooth and all their derivatives are bounded (C°°-bounded map); they represent 
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the no arbitrage conditions. Certainly ~ in view of existence and uniqueness - 
we have to assume global Lipschitz properties for the respective vector fields. 

2.1. Example. Our driving example will be a joint interest rate (IR) and foreign 
exchange (FX) market for p+l countries with different currencies and interest rate 
levels {p > 1). The risk factors of this market are the different log-exchange rates 
and forward interest rates. For the basic equations of this market we refer to [5]. 
Denote by Pt{x) and P^{x), a > 1 the time t prices of the domestic and foreign 
default-free, zero coupon bonds, where a — 1, . . . ,p denotes the different economies 
where exchange rates and zero coupon bonds are observed. The bonds are maturing 
dX T = t + X >Q. We define the domestic and foreign instantaneous forward rates 
as 

(2.5) r° = 91og P°, = d\og Pf. 

Using the Musiela parametrization and the standard Heath-Jarrow-Morton drift 
condition we can write down the domestic forward rate dynamics under the domes- 
tic martingale measure 

(2.6) dr? =(^r? -t- a(r°) . A. / a(r«) . K)dt + ^ a(rO) . A. dB^ 

i—l i—1 

Analogously under the foreign martingale measure we obtain a Heath-Jarrow- 
Morton equation for the foreign forward rate dynamics. If the dynamics of the 
spot log-exchange rate process with respect to the domestic martingale measure is 
as follows 

(2.7) dS^ = (r?(0) - r^rndt - (J^ ^-^)dt + J2 

i=l 2=1 

we obtain under the domestic martingale measure the foreign forward rate dynamics 
as 

(2.8) drf = (^rf + ^ a{r^) . A, J a{r^) . A.- 

d d 

i=i 1=1 

for a — 1, . . . ,p. Recall the no-arbitrage conditions, namely, the processes 

exp(- f r°(0)d<)F°(T), exp(- f r°(0)di) exp(S'")Pf(r) 
Jo Jq 

are (local) martingales on [0, T] for a = 1, . . . ,p and T > 0. 

This system of stochastic partial differential equations (SPDEs) is of the form 



(2.1 ), when we consider a vector of risk factors Y — (r", r^, . . . , r^*, 5^, . . . , S^) 



forming an element of a large Hilbert space 

(2.9) H := H° X ■■■ X HP xRx ■■■ xR. 

By abuse of notation we use the same letters Ai, . . . , A^ and a for the respective 
fields on the Hilbert space of risk factors, such that we obtain an equation of type 



(2.1 ). Notice that we can easily perform a change of measure in order to attain the 
physical measure. This corresponds then to introducing an appropriately modified 
drift, denoted by the additional term fi2iY). 
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3. Calibrating the diffusive SPDEs to the market 

In this section we propose a calibration to (diffusive parts of) the historical risk 
factor time series, which does not need neither optimization procedures nor the 
numerical evaluation of the model in question. Instead of using those standard 
approaches, which pose serious problems in high dimensions, we will rather con- 
sider a "layman's calibration" of the model, where the characteristics of the SPDE 
are deterministically and directly calculated from the observed time series. This 
corresponds to the "historical approach" , which is standard for several branches 
of scenario generation, however, we only apply the historical approach at an infin- 
itesimal level (thinking of a time-tick, usually from hours up to days in scenario 
generation, as the infinitesimal element): the global dynamics is then calculated by 
means of a stochastic differential equation whose local features are determined by 
the time series. 

We suppose that we are in the setting of Section [2] and we assume a time series. 



i.e. a sample of equation (2.11, on equidistant grid points of distance A, denoted by 
Yi, . . . , Yk- Our goal consists of estimating the volatility directions Ai, . . . , Ad out of 
the observations Yi, . . . , Yk in a simple way (we regard Yk as the last observation, 
which appeared one time tick A ago). We do really aim for an estimation in the 
proper sense, i.e., we would like to have a limit theorem for SPDEs that relates 
the number of observations K, the time-distance A of two contiguous observations, 
and the number of volatility directions d. We certainly do not want the dimension 
of the Hilbert space of risk factors to enter into the construction, since it might be 
infinite. 

For later use and for the sake of simplicity we state the stochastic partial differ- 



ential equation of type (2.1 1 calibrated to the given sample Yi, . . . , Yk 



(3.1) dxf^={^^^^\xl^^) + ^.[^\xf^))dt- 



K-l 

a(xf )) . {a{Y,)-\Y,+^ - y,)) dWl 



VA(X-l) 

where cr is a known, non- vanishing geometric function on the risk factors describing 
the local dynamics. We specify in the sequel the necessary conditions for an appro- 
priate limit theorem. We denote by x[^'^ the dependence of the solution process Xt 
on the number K of observations (which also determines the distance A between 
two observations) . The calibrated equation is defined on an independent stochastic 
basis with a, K — 1 dimensional Brownian motion {W^, . . . , W^~^). 

Notice here that ^[^"^ is specified by no-arbitrage conditions from the volatility 



structure, as given in the equations (2.7 1 and (2.8). Furthermore we have to assume 
the following technical assumption, which allows to construct strong solutions and 
not only mild (respectively weak) ones. For all necessary details see [B] and [T]. 

3.1. Assumption. The constant volatility directions Ai,...,Arf are elements of 
doTa(A°°) , which is the domain where all powers of A are defined (see for instance 
[1]. Furthermore a : dom(yi) dom(yi) is assumed to be a Lipschitz map. 



3.2. Theorem. Consider equation (2.1) where a is a given map and Ya G dom(yi) 
a given initial value, but the directions Xi,...,Xd are unknown. We assume As- 
sumption\3.1\ We collect a time series of observations Yi, . . . , Yk on an equidistant 



time grid of width A that cover an interval of length T — KA. Refining the obser- 
T_ 

K 



vations by making A = smaller and smaller leads to the following convergence 



statement: 

(3.2) lim - Yt 
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in distribution for any t > if Xq = Yq. The underlying limit result is the following 
Gaussian one, 

(3.3) lim f aixi^Y'dxl^^ - f aixfr\l^f\xi^^) + ^(XW)^. 

1 ^'^ 

(3.4) - lim V(a(y.)''(>;+i-y.))W^^ 

K^co ^A{K - 1) ^ 

(3.5) =/ <j{Yt)-'dYt- f a{Y,)-\^ii{Ys) + fi2{Ys))ds. 







Proof. Let H be the risk factor Hilbert space and X E H. Denote by A (K) A the 
bilinear operator defined by (A® A)(u, w) := (A, u) (A, lu), for {v,w) E HxH. Linear 
combinations of these (one-dimensional) bilinear operators are dense in the set of 
Hilbert-Schmidt operators (with respect to the Hilbert-Schmidt norm). Let M be 
a Hilbert space valued Ito process given by 



dMt = atdt + J2KdWi, 



i=l 

for a family a,b^ , . . . ,b'^ of square-integrable processes; then, the quadratic covari- 
ation 

J2 f{h\®b\)dt 

is approximated with respect to the Hilbert-Schmidt norm by 

K-l 
3=0 



see for instance [3] . The Gaussian random variable in equation ( 3.4 1 has a covariance 
matrix 

t ^"^ 

(3.6) A{K-1) ^^""^^'^ \Y,+i-Y,))®{a{Y^ \y,+, - Y^)) 

which precisely corresponds, up to the factor t/A, with the standard statistical 
estimator the returns (see Remark |4.1| for the use of the word return here) of 
the sample Yi, . . . ,Yk multiplied by (j{Yi) ^. On the other hand, this estimator 
converges, as K tends to infinity, to 
(3.7) 

K-l d 



^ ' i=l i=l 



This is due to the fact that the quadratic covariation of the stochastic process 

d 

(3.8) a{Yty'dYt = a{Yty\^lf\Yt) + fl2{Yt))dt + Y,^^dBl 

i=l 

is given by J2t=i{^i ® ^i) dt, whence ^ Ylt^ii^i ® ^i) is approximated with respect 
to the Hilbert-Schmidt norm by 

1 ^"^ 

A{K-l) ^^'"^^'^ \Y,+,-Yi))®{a{Y,) \y,+,-Y,)). 
Notice here that AK = T. This proves the equality of (3.3 1 and (3.5 1. 



SCENARIO GENERATION IN RISK MANAGEMENT 



7 



In order to prove the convergence result (3.2 1 we will have to apply the follow- 



ing limit theorem in distribution for a sequence of stochastic partial differential 
equations: let X^'^"' and Yt be the unique solutions of the equations 

(3.9) dX^ = (A^r^xf )) + ^i^\xi^^))dt + a{Xr^) . dBf\ X^ = ^o, 

(3.10) dYt = ifiiiVt) + fi2{yt))dt + <j{Yt) . dBt, Fo - ^o, 
where 

K-l 

(3.11) ) = V aiYr\Y.+i Wl 



d 

(3.12) Bt=J2>^^Bt 



for t > 0. Then 

(3.13) JCi'^' — > Yt 



in distribution as ii" ^ oo. In order to show (3.131, we start by noticing that, due 
to the previous consideration we have 

B(^) > B 

K^oo 



in distribution on the Hilbert space H. The limit theorem (3.131 can be concluded 
from the stability theorems in [5J and the uniqueness in law established in first, 
by uniqueness in law, the actual Brownian motion in use is not relevant. Hence, 
we can interpret the convergence of B^^^ as K ^ oo with respect to independent 
Brownian motions G\, each with a certain volatility tif^^ and associated to a certain 
eigendirection e|^\ that is, we can rewrite in law 

oo 

Br^Y.-f'ef'G\. 

1=1 

This yields, for each /, vector fields converging to a limit as if ^ oo, namely the 
eigensystem of Yl'i=i^i ® ^ii if o^i^ carefully enumerates the volatilities nf^^ in 
increasing order and keeps track of multiple eigenvalues. Second, if we reformu- 



late the stochastic partial differential equations (2.1) and (3.1 1 with respect to the 
sequence of independent Brownian motions G', we find ourselves with a sequence 
of vector fields converging to well-specified limits. Finally, the required statement 
follows from the stability result in Proposition 9.1 of [51. □ 

3.3. Remark. The methodology in the previous theorem is related to the one in [7] 
in the sense that we read from the quadratic variation process and its approxima- 
tions the relevant information for the a posteriori construction of a process which 
is close in law to the originally observed one. Indeed, in [7] the quadratic variation 
process and its approximations are - via harmonic analysis methods - investigated 
in order to reveal information on the non-linear volatility vector fields. In contrast, 
we do assume knowledge on the type of non-linearity in the equation and tackle 
the problem of determining the volatility directions Ai, . . . , A^. The convergence 
results based on Wiener's theorem, as stated in [f;, could be similarly applied. 

3.4. Remark. Notice that we do not need to calculate empirical covariance matrices 
of the time series in order to to capture the local correlation structure of all risk 



factors in the calibrated equation (3.1 1 



3.5. Remark. Up to some technical complications one could also consider d — oo 
here. 
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4. Rare events and Fine-tuning through random time changes 



Besides the calibration of the risk factor SPDEs proposed in equation (3.1 1 we 
have to take into account stochastic volatihty. We express stochastic volatiHty in our 
setting through the appearance of random times, which can be nicely interpreted 
as difference between physical time and trading time. As the spirit of this article 
is to embed a static sampling problem (in contrast to the distributional approach 
where only two random variables are described) into a dynamic problem defined by 
stochastic (partial) differential equations, we also embed the random time into a 
random time change, which can then a fortiori be observed through an stochastically 
changing volatility process. This approach is a continuous time version of GARCH- 
approaches, which are applied in filtered historical simulation in risk management 
(see [2]). Here, in contrast, we "calculate" stochastic volatility along the time series 
as realized volatility and rescale the time series by this realized volatility before 
we apply those rescaled ("filtered") returns to calibrate the underlying stochastic 



differential equations as proposed in equation (3.11 

In this section we therefore introduce the additional methodology, which allows 
us to include these changing volatility effects of the observed diffusion processes 
and which permits jumps. Both phenomena, stochastic volatility and jumps, are 
crucial to capture important stylized facts of financial time series. For instance, 
when we model rare (extreme) events we those events as additional jumps added 



to the diffusive setting of equation (2.11. Recall that these additional jumps also 
change the drift fii due to no arbitrage conditions. In this section we assume the 
existence of a basis {e^} of the respective Hilbert spaces of risk factors in order to 
speak about the risk factors , which denote the projection of Y onto the j-th 
basis element. 



4.1. Remark. We shall apply in the sequel the notion return for two slightly dif- 
ferent concepts. When considering a time series Yi, . . . , Yk, which is a collection 
of vectors, then we shall call any difference l^+i — Yi an (observed) return (of risk 
factors), where we should in fact speak of a vector of returns {Y^^^ of single 

risk factors Y^ . Running indices for risk factors are always denoted by j, the index 
along the time series is usually i. The length of the time series will be always be 
K. 



We define rare events using known time series of the risk factors and by the risk 
manager's very own intuition. Let Yi, . . . , Y^ denote a time series of risk factors. 
We form a time series of returns Yi+i — Yi for i = I, . . . , K — I from it and estimate 
the empirical (co-)variance a among its different components. Additionally we 
define local extractions of the time series of length L, i.e. Y^+i+i — Yk+i, for 
k = 1, . . . , L and for i — 0^ . . . , K — L — I. The empirical (co-) variance of this 
extraction is denoted by (t(*+^) and corresponds to the estimation of the volatility 
looking back L days, usually called (local) realized volatility. 

We now consider a diffusive equation of type (2.11 with jumps and random time 
change. More precisely, let = Aj^ds be a random time change for risk factor 
j, i.e. A' is a non-negative, locally integrable previsible stochastic process (with 
respect to the natural filtration of the Brownian and Poissonian componentes) , 
whence t l—^ Tj is an almost surely increasing random time. Then the equation to 
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be calibrated looks like 

d 

(4.1) dY^, - (Mi(y,J + ti2{Yr,))dri + (^(i^rj • X^YdBl, + dN^^ , 

' 1=1 

(4.2) dr/ = Aidt, 

(4.3) Yo e H, 

where N denotes a Poisson process with values in H and jump rate calibrated 
from the time series. The solution of this random SPDE can be written through 



the random time change by solving (2.1) up to the time r/ for the j-th risk factor. 
It is possible to estimate the volatility levels Al from the market up to an (a priori 
unknown) factor by observing a proper quadratic variation process between two 
jumps at time a and b, namely 

{[ {a{Y,)-'.dY,f)t= f AidsY, [\\ifds. 

As a second notation we introduce the ratios of local empirical volatility and 
today's local volatility, namely 



(K- 



which is defined for each risk factor j separately. We denote this quantity by R\ 
at time i in the time series. Notice that this quantity is well-defined, even under 
the regime of a random time change, since the (unknown) proportionality factor 
cancels out. 

4.2. Definition. Let Yi, . . . , Y/^' denote a time series of risk factors. We form a 
time series of returns Yi^i — Yi for i ~ I, . . . , K ^ 1 from it and denote by 

for risk factor j and i = L, . . . , K — 1 the rescaled returns of the time series. 

4.3. Remark. In order to make the choice of extreme events more stable along 
time, one can choose L for re-scaling the returns and for extracting extreme events 
independently. In the subsequent example this is the case. 



In order to calibrate the equation in the spirit of (3.1 1 we start from the set of 



rescaled returns and separate from it the set of extreme returns. 

4.4. Definition. A return at index L < i < K ~ I is said to be extreme at level rj 
with M violations if 

] _ v3 I 



YL. - Y: 



holds true for at least M different risk factors j. The set of all indices belonging to 
extreme events at level rj with M violations is denoted by E-n.M and is a subset of 
the set of all indices 3C = {I, . . . , K — I}. 

By means of the set we can now define a jump structure, where with a 

certain jump rate and jump measure jumps from £^^m are mixed (in an independent 
fashion) into the diffusion process calibrated with returns from the set 3C \ £ri,M- 

For the purpose of numerical efficiency of the simulation we shall most of the 
time neglect the fact that the random time change depends on each risk factor 
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. Actually the simulation time Tt is chosen with respect to the distribution of 
J^Asds, which is estimated from an averaged realized variance structure. 

5. Evaluating the stochastic differential equations 

In our case the parameters of the SPDE are specified by the the empirical cal- 
ibration method outlined in Section |3] and by an a priori choice of the geometric 
term cr. Therefore, in order to fully characterize the equation, the following points 
need to be taken into account: 

(1) Determine the sets £ri,M and 3C \ E.ri,M by fixing a level 77 at which events 
are considered as extreme and a number M of violations. 

(2) The number of involved Brownian motions is the cardinality of 3C \ E„m- 



(3) The volatility directions are given by rescaled returns using equation (3.1 1. 

(4) The semigroup A is the shift semigroup for term structure of the risk factors, 
otherwise the identity. 

(5) a is chosen a priori according to the (geometric) structure of the problem, 
i.e., some risk factors have to stay positive (like interest rates), others are 
free (hke log-prices). 



For the evaluation of the SPDE (3.11 we propose an explicit-implicit Euler scheme 
(see for instance [5 ), which has the following one-step Euler structure: 

• Prepare an initial risk factor Yq. 

• Trigger a jump (with a specified jump rate), or not. If a jump is triggered, 
then sample its size using the jump measure and add it to Yo- 

• Choose Brownian increments Bl with variance t along the time increment 
t. 

• Evaluate 

(5.1) StYo+^ii''\Yo)t + ^lf\Yo)t+ 

(5.2) + ' E ^iYo)'aiY.r^^^^^^-^Bl. 

• Terminate the Euler step. 

5.1. Remark. In our concrete implementations the choice of a depends on the 
respective risk factors. For instance, for interest rates we have chosen fields of 
linear operators depending via square-root like functions on the level of yield at 
certain maturities. 



6. Results from a concrete implementation 

In this last section we present results from a concrete implementation of a sce- 
nario generator for five foreign currencies (CHF, GBP, JPY, USD, ZAR) and six 
yield curves, i.e., five foreign yield curves and one domestic yield curve (Euro zone). 
For this purpose we use market data from March 2006 to March 2009, including in 
particular the highly volatile events of the last quarter of 2008. In total, this risk 
management problem consists of 82 risk factors, even though more points on the 
yield curve could have been evaluated. In a given day, the available historical time 
series has a length of K = 500 business days back in time, i.e., in March 2009 we 
look back up to March 2007, for instance. The implementation is based on Section 
|2] and the systems of equations presented therein. In the light of Sections |3] and 
Section [4] the realized variances (with L = 20) and the set of extreme events (for 
the calculation of extreme events we used L = 40) are calculated: the number of 
violations M — A and the level 77 = 4 yield a number of extreme events of 32 (for 
M — I, i.e. an extreme event occurs if we have one violation for one risk factor. 
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we obtain 162 extreme events). The jump rate is chosen to be 2 percent (yielding 
about 100 jump events in 5000 scenarios). The random time change has been cho- 
sen ad-hoc to be 0.9 days in 90 percent of the cases and 1.9 days in 10 percent of 
the cases for the sake of simphcity. A precise cahbration might yield even better 
results. 

The following pictures show several single distinguished historical and gener- 
ated distributions of risk factors. The historical distributions are rescaled by R 
as described in Definition |4.2[ The last figure shows the backtesting results for a 
randomly chosen roll-over portfolio of linear FX- and IR-instruments, i.e., the times 
to maturity of instruments in the portfolio are kept constant along the backtesting. 
The distributional plots show on the left hand side the historical distribution of 
log-returns for one day and on the right hand side the simulated distributions of 
log-returns for one day. 

It is apparent that our simulated scenarios resemble the historical distributions 
by smoothening them in a consistent way along all risk factors and by preserv- 
ing extreme risks. We obtain backtesting results of four VAR-violations and two 
portfolio-loss-values beyond the level of expected shortfall in 250 days at a 99% 
confidence level. 

In the first set, the smoothening is the main visible effect in the simulation 
results. Also other stylized facts of the distribution like skewness and kurtosis are 
preserved. 
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In the second set of plots, mass of the historical distributions is transported from 
the center to middle-size events. 
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In the third set of plots the mass transport effect is seen in an even more extreme 
way, since for ZAR a lot of returns have historically been zero. In this sense the 
constructed scenario generator also has a "repairing" effect on data. 
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Finally, we present backtesting results of a two roll-over portfolios during 250 busi- 
ness days at the 99% confidence level. Returns are in red, VaR in blue and the 
value of expected shortfall is in green. One can clearly see the increase of the level 
of VAR (and expected shortfall) during the events of the last quarter 2008. We 
have omitted the scaling on the y-axis as violations are only a level of relative size. 
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